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Abstract 

We justify some characterizations of the ground states of spin-1 Bose-Einstein con- 
densates exhibited from numerical simulations. For ferromagnetic systems, we show the 
' (— I '. validity of the single-mode approximation (SMA). For an antiferromagnetic system with 

Oh! nonzero magnetization, we prove the vanishing of the mp = component. In the end of 

(~| ' the paper some remaining degenerate situations are also discussed. The proofs of the main 

"j^ ■ results are all based on a simple observation, that a redistribution of masses among different 

components will reduce the kinetic energy. 

> : 1 Introduction 

(N 
m 

00 ! At ultra low temperature, massive bosons could occupy the same lowest-energy state and 
[ form the so-called Bose-Einstein condensates (BECs). This phenomenon was predicted 
§ ! by Bose and Einstein in 1925, and was first realized on several alkali atomic gases in 
1995 by laser cooling technique [HIEIISI. In early experiments, the atoms were confined 
in magnetic traps. In this situation the spin degrees of freedom are frozen. Through 
the mean-field approximation the system is then described by a scalar wave function, 
which satisfies the Gross-Pitaevskii (GP) equation [l8l[T3l|2ll]. In contrast, in an optically 
trapped atomic BEC all hyperfine spin states can be active simultaneously, and a spin-F 
BEC is then described by a vector wave function \E' = {ipF,'ipp_i, ■ ■ ■ , ip^p)'^, where the 
j-th component corresponds to the m^? = j hyperfine state [|24l |251 [191 IH [T2]] . The theory 
of such spinor BEC was first developed independently by several groups |l20l [IH [T6l . 
After these early studies, spinor BEC has become an area of great research interest. 
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1.1 Mathematical model for spin-1 BEC 



For a spin-1 BEC, the vector wave function 
GP equation: 



6E 



V'Oj'^-i)"^ satisfies a generafized 



(1-1) 



where the Hamihonian is given by 



3 + £2|*|4 + ££|$-s$|2 



>dx. 



Here D is a domain in M.'^, h is the reduced Planck constant, is the atomic mass, \^ is a 
locally bounded real- valued function representing the trap potential, \&* is the Hermitian 
of ^, and S = {Sx-,Sy,Sz)'^?' the triple of spin-1 Pauli matrices: 



Sx — ^ 
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So ^*S'^ denotes the vector (^*5'a;^, ^*S'y^, ^*S'2^). Also note that |^| denotes the 
Euclidean length l^jP)^^^' ^^'^ similiarly for |VV^j| and |\E'*S'\E'|. The parameters Cn 
and Cs Sire real constants given by 



3m„ 



(ao + 2a2), Cs 



47r/i2 
3ma 



{-ao + 02), 



where ao and 02 are respectively the s-wave scattering lengths for scattering channels 
of total hyperfine spin zero and spin two. The parameter characterizes the spin- 
independent interaction, and the parameter characterizes the spin-exchange interac- 
tion. For Cn < (resp. Cn > 0), the spin-independent interaction is attractive (resp. 
repulsive). For Cg < (resp. Cg > 0), the spin-exchange interaction is ferromagnetic 
(resp. antiferromagnetic). Typical examples of ferromagnetic and antiferromagnetic 
systems are ^''Rb and ^"^Na condensates. 

The generalized GP equation (11.11) implies two conserved quantities: 



(CI) 
(C2) 



D 



D 



(|V^l|2-|^„l|2) =M, 



where N is the total number of atoms and M is the total magnetization. For the system 
to be nontrivial, we assume > 0. We also assume \M\ < N (note that obviously 
1^1 < N), for if \M\ = N the system reduces to a single component BEC, which is a 
trivial case for all considerations in this work. Now we say ^ is a ground state if it is a 
minimizer of E under the above two constraints. 
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1.2 Innovation and organization 



In researches concerning ground states of spin-1 BEC, the following ansatz was often 
adopted: 



where Cj are constants and ip is a. function independent of j. This is called the single- 
mode approximation (SMA) in the physics literature [US [IH |2l [Bl |22l HQl . It has been 
found f2E\ from numerical simulations that ground states obey the SMA exactly for fer- 
romagnetic systems (and does not in general for antiferromagnetic ones), and hence can 
effectively be characterized as one-component systems. The first goal of this paper is 
to analytically confirm this observation. On the other hand, for antiferromagnetic sys- 
tems, we will show that ipo = when M 7^ 0, another well-known phenomenon from 
numerical simulations |j3ll21 not being rigorously proved before. For the degenerate case 
M = 0, however, the SMA is again valid while ground states are not unique, and ipo does 
not necessarily vanish. It's interesting that although the two phenomena (SMA and van- 
ishing of tpo) look quite irrelevant to each other, they can be proved by the same simple 
principle, that a redistribution of masses between different components will decrease the 
kinetic energy. 

The paper is organized as follows. Section 2 is the preliminary, where we reformulate 
the mathematical model more precisely, and then provide a result of maximum principle 
which is crucial in justifying the expected characterizations. In Section 2.2 the idea of 
mass redistribution is introduced. Sections 3 and 4 treat respectively the ferromagnetic 
and antiferromagnetic systems. 

2 Preliminary 

For notational simplicity, let's redefine 



This causes no loss of generality for the phenomena we are going to investigate. The 
admissible class is 



where L?'{D,Vdx) consists of all functions / such that j^V\f\^ < 00. Let u denotes 
Mo, We also define 

^ = {u G C I > for each j] ; 

^1 = {u G ^ I u = (71/, 70/, 7-i/) for some constants 7^ and some function /} ; 
^2 = {u G ^ I no = 0} . 



cjip for each j, 
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Let's also use 7 to denote (71, 70, 7-1), so (71/, 70/, 7-1/) can be abbreviated as 7/. 

In Section 2.1, we introduce a common reduction which shows that to study ground 
states we can simply consider A instead of C. Indeed, A consists just the amplitudes of 
elements in C. And Ai (resp. A2) corresponds to the set of all elements obeying the 
SMA (resp. with vanishing zeroth components). For the moment, we do not consider 
any boundary condition for simplicity. See the remark in the end of Section 3. 



2.1 Reduction from C to A 



Given \I> g C. Let uje"^^^ be the polar form of ipj. It's easy to check that if \E' is a ground 
state, that is minimizes E over C, then the Oj's are constants satisfying 



cos (^1 - 200 + 9^1) = ±1 for cs ^ 0, 



and 



E[^] 



D 



J2 1^".^ 



Vlup + Cnlul'^ + Cs 



(2.1) 



(2.2) 



where the plus-minus sign ± corresponds to Cg ^ 0. Let's now define E : ^ — M, E[u] 
is given by the right-hand side of (12.21) . What we claimed is if is a ground state, then 
E[if] = K[{\ipi\, lipol, \ip^i\)]. Conversely, if any u e A satisfies 



El 



minE[v] 
veA 



the vector \1/ defined by ipj 



u 



,e*^i is a ground state as long as the ^/s are constants 



satisfying (12.11) . Thus, studying ground states of E is equivalent to studying minimizers 
of E. Without loss of generality, we will henceforth consider E instead of the original 

E. 

For convenience let's use H to denote the integrand of E, i.e. E[u] = H{\i). We 
also write H = Hi + H2, where 



Hi{u) = ^ iVwjf + Cs 2ul{ui ± u-if + (uj 



u_ 



i72(u) = V|u|2 + Cn|u|l 



This splitting of H has no physical meaning but only for convenience of later discussion. 

We shall denote the set of all minimizers of E over Ahy Q. The Euler-Lagrange 
system for u g ^ is given by the following coupled Gross-Pitaevskii equations: 



(/i + X)ui = Cui + 2cs [uiiui ± + ui{u\ - 1)] 

HUq = Cuq + 2CsUq{ui ± 

(/i - A)m_i = + 2cs [ul{u^i ± Ml) + u^i{u^i - -uf)] 



(2.3) 
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where C. = —A + V + 2cn\u\'^, A and /i are the Lagrange muMpliers. We remark that in 
this paper we do not involve ourselves in the problem of existence. To best illustrate the 
simplicity of our method, we just assume there is a ground state, (see [[TSlOEl] for related 
concerns of existence problem). Also note that u g ^ is continuously differentiable by 
standard regularity theorem (see e.g. QlTJ 10.2]). 

The following lemma will be crucial in our characterizations of ground states. 

Lemma 2.1. IfueQ, then for each j, either uj = or uj > on all of D. 

Proof. For an arbitrary compact K c D, by subtracting respectively QjUj, j = 1, 0, -1, 
from the three equations in (12.31) with large enough constants Qj, and using the assump- 
tion Uj > 0, it's easy to verify that each Uj satisfies 

Auj + hjUj < 

for some hj < on . Thus either uj > or = on by strong maximum principle. 
Since K c D is arbitrary, the assertion of the lemma holds. □ 



2.2 A kinetic-energy-reduced redistribution 

Consider an n-tuple of nonnegative functions f = (/i, /2, fn) e {H^{D)Y. Let g = |f |. 
It's well-known that \Vg\^ < l^/fcP- In f'^ct, direct computation gives 

, -2 Sl/jWfc - /fcV/jf onwhere^>0 
5]|V/fcP-|V(7|2= <{ 3 (2.4) 

) on where g = 

In particular, \Vg\^ = Y.k \^fk? if and only if fjVfk - fk^fj = for j ^ k. 

The property above has a simple while interesting generalization, when /| {k = 
1, 2, n) do not sum to a single g'^, but instead are redistributed into multiple parts. To 
be precise, we give the following definition. 

Definition 2.1. Let f be as above, and let g = {gi,g2,..., gm) be an m-tuple of nonnegative 
functions. We say g is a mass redistribution of f if 

n 

fc=l 

for each £ = 1, ...,m, where a^k are nonnegative constants satisfying YlT=i = 1, for 
each = 1, ...,n. 



Note that = |f | is the only mass redistribution of f with m = 1. In general we have 
the following proposition. 
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Proposition 2.2. For any mass redistribution gofi as in Definition 12. II we have 



(1) |g| = |f|; 

(2) J2T^i l^ail"^ < YTk^i I V/fcP- Equality holds if and only iffjVfk - fk^fj = Ofor 
each j ^ k with aijaik ^ Ofor at least one I. 

Proof. The first assertion follows directly from the definition of redistribution. For (2), 
from (12.41) we have 

\ ^ aijaiklfj^fk - fk^fj? on where > 

I on where ge = 0, 

and the assertion is obtained by summing over all i. □ 

In this work, we will consider mass redistributions of u g ^. Since the square of 
the amplitude of a wave function represents its mass distribution, we use the adjec- 
tive "mass" to stress that it's a redistribution of the square^. To save notation, in the 
following we shall omit it and simply say "redistribution". Note that if u g ^ and 
V = {vi,VQ,v^i) is a redistribution of u, then from (1) of Proposition 12. 2[ v satisfies 
(CI) automatically and H2iv) = H2{u). These facts together with (2) of the proposition 
allow us to give unified and simple justifications of the two properties mentioned in the 
introduction. 



3 Ferromagnetic systems 

In this section we assume Cg < 0, and the goal is to prove the validity of SMA. That is we 
want to show Q c Ai. The idea is to find, for u g ^, a redistribution of u in Ai that has 
no larger energy than u, and then try to conclude that u must itself be the redistributed 
element provided u e Q. 

Now given any u e A. It's easy to see that a redistribution of u in Ai can be expressed 
as 7|u|, where 7 = (71, 70, 71) is any triple of nonnegative constants satisfying 

7? - 72 ^ = M/N. 
Let r denote the set containing all such 7: 

r := {7 G I 7j- > for each j, 7 satisfies (El])} . 



'And one can naturally generalize the idea to p-th power redistribution, which may be useful in studying 
systems with p-Laplacian terms. 
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Then 7/2(7! u|) = ^2(u) for each 7 g F. On the other hand, 

i^l(7|u|) = |V|u||2 + c,P(7)|u|^ 

where 

M2 



P(7) = 27o"(7i + 7_i)^ + 



Ar2- 
Since Cg < 0, 

min (c5P(7)|u|^) = c, I maxP(7) ) |u|'*, 
and it's easy to check that 

maxP(7) = P(7*) = l, 
where the maximizer 7* = (71 , 7o, 7* 1) is unique and is given by 



On the other hand, we akeady know | V|u| p < I Vwjp. Thus we have proved 

H{j*\u\) = |V|u||2 + cs\u\^ < H{u) (3.2) 

for any u e A. We can now prove the first characterization of ground states by examining 
the condition for equality of (13.21) . 

Theorem 3.1. Assume Cs<O.IfueQ, then u = 7*|u|. 

Proof. 

H{u) - H{-f*\u\) = Hi{u) - Hi{-f*\u\) 





|V|u|| 




|V|u|| 



by direct calculation. If u g ^, we must have E[u] = E[7*|u|], which implies H{u) = 
if (7*|u|), and hence 

UjVuk — UkVuj = for i k] (3.3) 
ul - 2uiu-i = 0. (3.4) 

Since we assume the total number of atoms > 0, from Lemma IXTI at least one Uj is 
strictly positive in D. Assume mi > on D. Then from (13.31 ) we have 



Since D is connected, (13.51) implies no and n_i are both constant multiples of ui. This 



shows u e Ai, and (ii) follows either by (13.41) or by the fact that 7* is the unique maxi- 
mizer of P over T. The case mq > and m_i > can be proved similarly. □ 

The above theorem implies that searching for ground states of ferromagnetic spin-1 
BEC can be reduced to a "one-component" minimization problem. Precisely, let 

= {|u| I u G ^} = {/ G H\D) n L^{D) n L^{D, Vdx) \ Jd f = ^} ^ (3-6) 
and define : ^ R, 



f {|V/|2 + V/2 + (cn + c,)/4}. 

JD 



W[f] 

Then E[Yf] = E'[f] for f eA'. Also let 

g'=lfeA' W[f] = mmW[g]]. 

Then if u g ^, u = 7*|u|, and hence 

E^[|u|]=E[7lu|]<E[7V]=E^'[/] 
for every / g A^. Thus |u| g Q^. Conversely if / G then 

E[7V] = mf] < mM = E[7iu|] < E[u] 

for every u e A, and hence 7*/ G Q. We thus obtain the following characterization of Q. 
Corollary 3.2. g = {7*/ | / g G'}. 

Remark. We can add more assumptions in the definition of A. The only thing we need 
to take care is that we need 7*|u| g A whenever u g ^, so that E[u] < E[7*|u|] is not 
violated. In particular, in case that a homogeneous boundary condition (e.g. homoge- 
neous Dirichlet or Neumann boundary condition) is considered, the induced boundary 
condition for 7*|u| is also homogeneous of the same kind, and Theorem 13. II (and hence 
Corollary 13.21 ) remains valid. 



4 Antiferromagnetic systems and some degenerate cases 

The main focus of this section is the phenomenon uq = 0. After justifying it in Section 
4.1, some degenerate situations are also discussed in Section 4.2. 
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4.1 Justification of tlie vanishing phenomenon 

Assume > in this subsection. We want to show that any ground state must have 
a vanishing zeroth component, and hence is a two-component BEC. Similar to the ap- 
proach in the previous section, we want to find an appropriate redistribution u g ^2 of 
u G ^ so that E[u] < E[u]. Now, not as before, the assumption u e A2 doesn't give rise 
to a definite candidate of u. In view that such u satisfies |u| = |u| and hence (CI), as a 
guess, we try imposing the additional assumption that u also satisfies 



~2 ~2 2 2 

Ml — U_-t = U^ — U_ 



1) 



SO that (C2) is also satisfied by u automatically. This results in only one possibility, that 
is 

«i = y«2 + !| for J = 1,-1. (4.1) 

For any u g ^, we then let u g ^2 be its redistribution defined by (14.11) . It's fortunate 
that it works. In fact, 

H{vi) - H{u) = ( ^ iVujf - ^ iVwjP) + 2csul{ui - u^if > 0, (4.2) 
j j 

and we have the following analogue of Theorem 13. 1[ 

Theorem 4.1. Assume Cs > and M ^ 0, then u g ^ implies u = u. 



Proof. Assume u g ^. From (lO) . ]E[u] = ]E[u], and 

ul{ui - u^if = 0. 

By Lemma [Xn we have either mq = or mi = m i. However since we assume M ^ 0, 
we cannot have ui = and the assertion follows. □ 

4.2 Some degenerate situations 

The requirement M 7^ in Theorem 14.11 is necessary. In fact, for M = 0, SMA is 
again valid while ground states are not unique, and mq = is not necessarily the case. 
Precisely, consider the minimization problem (recall that A'^ is defined by (13.61) ) 

min ! {|V/|2 + r/2 + c./4}, (4.3) 

and we have the following characterization. 
Theorem 4.2. 7/"cs > 0, M = or = 0, then 

^ = I (t, - 2t2,t)/ < t < 1/^2, / is a solution o/dS])} • 
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Proof. Note that since M = 0, from (13.11) . 7 G T implies 

7 = (t, - 2^2, for some t g [O, 1/^2] . 
Now it's easy to see that for any ue A and 7 G F we have 

iJ(7|u|) = |V|u||2 + l^|u|2 + Cn|u|^ 

which satisfies if(7|u|) < H{\i) by Proposition 12.21 and the remaining of the proof is the 
same as in Section 3. □ 

In contrast to the above theorem, the following corollary of Theorem 14. 1 1 shows that 
SMA is almost never the case when M 7^ 0. 

Corollary 4.3. Assume > and M 7^ 0, then u g ^ fi ^1 implies ui, u-i and V are 
constants. 

Proof. By Theorem 14.11 the Euler-Lagrange system (12.31 ) is reduced to the following 
two-component system: 

+ X)ui = Cui + 2csUi{u\ - u\) ^^^^ 
{fi — X)u-i = Cu^i + 2csU-i{u^_i — Ui), 

where C = -A + V + 2cn{u\ + u\). 

Recall that we assume -A^ < M < N, thus, for j = 1,-1, Uj > on D. So u g ^1 
implies m 1 = kui for some constant k > 0. Also note that k ^ 1 since M 7^ 0. The 
system (|4.4|) then gives the following two equations for m: 

(yu + \)ui = -Ami + Vui + 2cn(l + k'^)uI + 2cs(l - k'^)uI; (4.5) 
(/i - X)ui = -Ami + Vui + 2cn(l + k^)uI + 2cs{k'^ - l)ul. (4.6) 

Now (14.51) minus (14.61) gives Xui = 2cs(l - k^)uI. Since wi > on D, we get 



/ A 

Ml = 



. 2c,(l-ft:^) 

In particular ui and m_i = nui are constants. Hence Ami = 0, and then (|4.5I) plus (|4.6I) 
gives 

/iMi = VUl + 2Cn(l + 

from which we get 



= /X - 2Cn(l + K^)Mi = /i A, 

Cell K 



which is also a constant. □ 
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